A General Theorem Relating the Bulk Topological Number to Edge States 
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We prove a general theorem on the relation between the bulk topological quantum number and 
the edge states in two dimensional insulators. It is shown that whenever there is a topological order 
in bulk, characterized by a non-vanishing Chern number, even if it is defined for a non-conserved 
quantity such as spin in the case of the spin Hall effect, one can always infer the existence of gapless 
edge states under certain twisted boundary conditions that allow tunneling between edges. This 
relation is robust against disorder and interactions, and it provides a unified topological classification 
of both the quantum (charge) Hall effect and the quantum spin Hall effect. In addition, it reconciles 
the apparent conflict between the stability of bulk topological order and the instability of gapless 
edge states in systems with open boundaries (as known happening in the spin Hall case). The 
consequences of time reversal invariance for bulk topological order and edge state dynamics are 
further studied in the present framework. 

PACS numbers: 73.43.-1,71. 10.-w,71.10.Pm,72.25.Dc 
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I. INTRODUCTION 

Since the discovery of the quantized Hall effect two 
decades ago 1 -, topological quantum number has been well 
accepted as an important characterization of quantum 
many-body systems. The quantized Hall conductance 
in the integer quantum Hall effect is first identified, for 
non-interacting electrons on a lattice, as a first Chern 
number over the magnetic Brillouin zone by Thouless et 
alrk, which is also known as the TKNN number. How- 
ever, this characterization does not apply to systems with 
impurity and/or interactions, since the Brillouin zone is 
not defined for these systems. The study on momentum 
space topology 3 can be viewed as the developments along 
this direction. It indeed goes beyond the non-interacting 
fermion system, by considering the poles of the single par- 
ticle Green function, but it pre-requires the Fermi liquid 
property of the system. To define a topological char- 
acterization for general many-body systems, Niu et alX 
introduced the technique of twisted boundary conditions, 
and showed that the first Chern number can be gener- 
ally defined over the parameter space of twisted phases 
for any two dimensional insulator. 

One of the important observable consequences of the 
bulk topological invariance in quantum Hall systems is 
the presence of gapless (chiral) edge states on the bound- 
ary of the system. Just like the existence of gapless Gold- 
stone modes is a characteristic of the spontaneous break- 
ing of a continuous symmetry, the existence of gapless 
edge excitations can be considered as a manifestation 
of bulk "topological order", defined through the topo- 
logical quantum numbers. According to Laughlin and 
Halperin's gauge argument* 5 ^, n branches of gapless edge 
states must exist for a system with Hall conductance 
o~h = ne 2 /h, since there must be n electrons transferred 
from one edge to the other when a unit flux 2tt is adiabat- 



ically threaded through a cylindrical system. The phys- 
ical intuition of this argument inspired later Hatsuga: 7,8 
to demonstrate an explicit and rigorous relation between 
the bulk TKNN number and the dynamics of edge states. 
He has succeeded in finding a topological characterization 
of edge states and then proved that the edge topologi- 
cal number is equal to the bulk TKNN number. How- 
ever, this relation cannot be generalized to more general 
many-body systems, where the TKNN number has to 
be replaced by a more general Chern number as defined 
in RefA Therefore, it is important to establish a gen- 
eral theorem, similar to the "Goldstone theorem" in the 
case spontaneous symmetry breaking, that bridges the 
bulk topological quantum number generally defined by 
the twisted boundary condition, and edge properties for 
a generic many-body system. To provide such a theorem 
is the main goal of the present paper. 

This problem is not merely of academic interest, in 
view of the recent proposals of the intrinsic spin Hall 
effect 9 ' 10 : 11 ' 12 and the quantum spin Hall effect 13 : 14 ' 15 . 
Two distinctly new issues arise in the topological char- 
acterization of the quantized spin Hall effect, compared 
to that of the quantized charge Hall effect, calling for 
a clearer understanding of the general relation between 
bulk topology and edge dynamics. Firstly, spin is in gen- 
eral not conserved due to spin-orbit coupling and, there- 
fore, the Laughlin-Halperin gauge argument does not 
apply straightforwardly. However, by generalizing the 
twisted boundary conditions in Ref.— to the spin chan- 
nel, a Chern number can still be definecUSiiLi^. Thus 
the physical meaning of such bulk Chern number and its 
relation to the edge states in the quantized spin Hall case 
becomes less clear compared with the quantum Hall case. 
Secondly, the spin Hall systems of present interest respect 
time-reversal invariance, which leads to additional con- 
straints on the possible Hamiltonians and on the edge 



dynamics. In Refii£, Kane and Mele proposed a "Z2 
topological order" for T-invariant systems, in which the 
T-invariant systems are classified according to whether 
there are an even or odd number of Kramers pairs of edge 
states on each edge. Such a Z2 characterization is pro- 
tected by Kramers degeneracy and seems to be consistent 
with recent numerical resultsA£*A&. However, it remains 
unclear whether the Z^ characterization can be general- 
ized to many-body systems with impurities and/or inter- 
actions. Careful study on edge dynamica2°*2I shows that 
the edge states can become gapful with sizable strength 
of interactions, even for the case of an odd number of 
Kramers pairs. This suggests that the Z% classification 
may not be a topological order protected by the bulk 
gap alone, but depends on more subtle behavior of the 
system. 

To clarify these issues, in this paper we propose a new 
framework to relate the bulk topological number to edge 
dynamics, by introducing generalized twisted boundary 
conditions. In the usual definition of twisted boundary 
conditions, the system is defined on a torus (namely a 
rectangle with side length L with the opposite boundaries 
in x- and y-directions identified, respectively); a particle 
gains an extra phase e l9x (or e v ) whenever it goes across 
the boundary in x (or y-) direction. (The phases (9 X , 9 y ) 
may be associated with the internal degrees of freedom 
in the case of the spin Hall effect.) Assuming that the 
many-body ground state of the system is separated from 
the excited states by a finite gap for all values of the 
twisted phases (9 x ,9 y ), the first Chern number is then 
defined by the total flux of the Berry-phase gauge field 
associated with the ground state over the (6 X , 6 y )-space, 
whose topology is a torus. To establish a bulk-edge con- 
nection, one may relate the toroidal system to a cylindri- 
cal system with open boundaries in the way as suggested 
in Reff2£. The simplest thought is of course to cut the 
torus in real space along a circle, say located at x — L, re- 
sulting in a cylinder with two open boundaries at x = 
and x = L. Then a particle can no longer go across 
the boundary at x = or x = L. Conversely, one may 
view the toroidal system as obtained by gluing together 
the two edges of a cylindrical system. In quantum theory 
the gluing can be implemented by allowing a tunable tun- 
nelling between the two edges2£. This leads to the idea 
of changing the twisted phase e 10 * of the toroidal system 
into a complex number re l6m , with < r < 1. Namely 
we will allow the hopping amplitude of a particle across 
the x-boundary to be reduced by a factor r. When r = 0, 
no particle can go across the x-boundary, and the spa- 
tial geometry of the system becomes a cylinder. With 
this generalized twisted boundary condition, the param- 
eter space is now enlarged to three variables (r,9 x ,9 y ). 
Since at r — 0, different values of 9 X do not make any 
difference, the geometry of the parameter space is actu- 
ally D 2 x U(l), with (r,9 x ) the polar coordinates of the 
disk D 2 . In other words, the topology of the parameter 
3-space is now that of a solid torus, the interior of a torus 
embedded in three dimensions: The toroidal surface at 




FIG. 1: Schematic picture of the parameter space torus and 
the definition of coordinates (r, 9 x ,0 y ). A monopole singular- 
ity is drawn in the solid torus. 



r = 1 represents the original toroidal system, while the 
central loop at r — represents the associated cylindrical 
system. In this way, the solid torus implements a contin- 
uous interpolation between the toroidal and cylindrical 
systems. A schematic picture of the parameter space is 
shown in Fig. ^ In this framework, the bulk Chern 
number corresponds to the total magnetic flux through 
the toroidal surface at r = 1, so inside the solid torus 
there must be magnetic monopole sources, i.e. singular- 
ities where degeneracy occurs due to the emergence of 
gapless edge states and breaks the adiabatic evolution 
and makes the Berry phase ill-defined. Consequently, 
the relation between bulk topology and edge dynamics is 
simply dictated by the (magnetic) Gauss theorem for the 
solid (r,9 X7 9 y ) torus. 

Once such a novel connection is established, more 
complete understanding on the quantum Hall (QH) and 
quantum spin Hall (QSH) systems can be achieved. It 
will be demonstrated that whenever the bulk topology 
is nontrivial, there must be gapless edge states in some 
twisted-boundary system with (r,9 x )(r < 1), but the 
gaplessness of edge states does not have to occur in the 
open-boundary system with r — 0. Only when the corre- 
sponding physical quantity in the twisted boundary con- 
ditions (charge for QHE and spin for QSHE) is conserved, 
can the gapless edge states be sure to happen in the open- 
boundary system. The "shift" of gapless edge states to 
generalized twisted boundary conditions provides a nat- 
ural rationale to account for the gap opening in the edge 
channels in the spin Hall systems22i2i. Upon consider- 
ing the time reversal invariance, some constraints on the 
winding number distribution is obtained and the Z<i clas- 
sification is shown to be related to the behavior of low 
energy excitations under time reversal, not merely to the 
bulk topology. 

The rest of the paper is organized as follows: The 
proof of the general theorem on the relation between bulk 
topology and edge dynamics is given in Sec. II. Several 



examples are given in Sec. Ill, including the QHE and 
QSHE cases. The role of time reversal invariance is dis- 
cussed and the Z2 classification is analyzed in Sec. IV. 
Finally, Sec. V is devoted to summary and discussions. 



II. GENERAL RELATION BETWEEN BULK 
TOPOLOGY AND EDGE STATES 

A. Generalized twisted boundary conditions and 
the bulk Chern number 

There are two equivalent ways to define the twisted 
boundary conditions for a toroidal system. One is to 
impose constraints on the many body wave function of 
the form $(ri,r 2 , ..) = e i6, $(ri +Lx,r 2 , ..). The other is 
to add parameters into the Hamiltonian. Since we would 
like to study the Berry phase gauge field under adiabatic 
evolution, it is more convenient to choose the latter. For 
convenience, here we consider the many-body system on 
a lattice, with a Hamiltonian of the form 



Ho 



(ij) 



Hu 



(1) 



in which Ci a is the annihilation operator of boson or 
fermion at the site i, with a = 1, • • • , s labelling internal 
degrees of freedom. The first term is the nearest neigh- 
bor hopping term. The second term H lnt may contain 
all possible disorder or interaction terms that conserve 
the local particle number: [H- lntl n.i\ — 0,V«. The sites 
i are defined on a 2-d lattice with size L with periodic 
boundary conditions or, in other words, a lattice on a 2- 
torus T 2 . Such a Hamiltonian (Q describes a wide class 
of interacting many-particle systems with nearest neigh- 
bor hoppings. Our treatment below can be generalized 
to more complex Hamiltonians as discussed later. 

We introduce on T 2 two boundary lines L x , L y in the 
x and y direction, respectively, as shown in Fig. [3 The 
twisted Hamiltonian is defined by changing the hopping 
matrix elements across the boundary lines L x , L y only in 
the following way: 
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with the parameters r, 9 X and 9 y being real. Here the 
hopping across boundary lines is referred to that from 
the left to the right or from the down to the up, as shown 
in Fig. [3 For the hopping across the boundaries in the 
opposite direction, we have to take complex conjugate 
to keep the Hamiltonian hermitian. All other hopping 
that do not cross L x or L y remain unchanged. If r = 1 
and r^ and T y are the unit matrix in internal indices, 
Eq. J2J describes the usual twisted boundary conditions 
introduced in Ref^, with 8 X and 9 y the twist phases due 
to flux threading. 

The above boundary conditions @ generalize the 
usual ones in two important aspects. First we have intro- 



duced generators T x and T y that act on the internal in- 
dices. In the spin Hall effect, they can be used to describe 
spin-dependent twisted phases^. We need to normalize 
the generators so that all their eigenvalues are integers 
so as to make the phase factor e "M^M 27r-periodic 
in 9 Xl 8 y . The more important feature of our boundary 
conditions @ is that we have introduced a real factor r 
that modifies the amplitude of the hopping across the x- 
boundary. With r = 1 one recovers the original toroidal 
system, while with r = the particles can not go across 
or beyond the x-boundaries so that the system becomes 
a cylindrical one, having open boundaries. Obviously for 
the values < r < 1, we allow a tunnelling between 
the two edges of the cylindrical system. In this way, our 
twisted Hamiltonian is parameterized by three real pa- 
rameters r, 9 X , 9 y and two hermitian generators T x ,T y : 



H = H, 



r„r, 



(r,( 



(3) 



From Eq. J5J it is clear that if r — all values of 9 X are 
equivalent to each other, so the topology of the (r, 9 X , 9 y )- 
space, with < r < 1 and < 9 X , 9 y < 2ir, is a three di- 
mensional solid torus. Our twisted Hamiltonian describes 
a continuous interpolation between a cylindrical system 
with r = and a toroidal system with r = 1. The imple- 
mentation of this interpolation is, as we will see below, 
crucial to establishing a general theorem that relates the 
bulk topological quantum number of a toroidal system 
to the edge dynamics of the associated open-boundary 
system. 

Before proceeding we make the remark that more com- 
plex hopping terms, like next nearest neighbor hopping 
etc or even two-particle hopping cJ a cLcj- 7 Cjjt^a„^, can 
be included in the Hamiltonian, as long as all the hop- 
pings are of a finite range. The twisted boundary con- 
ditions are defined by a straightforward generalization of 
Eq. (01 for the hoppings across the boundaries. All the 
discussions below will remain essentially unchanged. 

With the parameterized Hamiltonian, the Berry phase 
gauge field can be defined at any point (r, 9 x ,9 y ), where 
the Hamiltonian has a unique ground state with an exci- 
tation gap. As usual the Berry gauge potential is given 
by 

Ap(r, 9 x ,9 y ) = -i (G(r, 9 x ,9 y )\ J- |G(r, 9 x ,9 y )) , (4) 

with fi — x,y. The system acquires a Berry phase under 
an adiabatic evolution along loop C as the loop integral 
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(5) 



For a toroidal system in an insulating phase, we take 
r = 1 and assume that the Hamiltonian -ffr^r (Ij $x> @y) 
has a non-degenerate ground state for all 8 x ,9 y . The 
bulk first Chern number is defined as the total flux of 
the gauge field through the (9 x ,9 y )-torus with r = 1: 
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FIG. 2: Schematic picture of twisted boundary conditions J2J. 
The a: and j/ boundaries L^ and L y , respectively, are drawn 
in red dash lines. Each hopping across them in the direc- 
tion indicated by the arrows acquires an additional matrix 
factor re and e y y , respectively. The hopping in the 

opposite direction acquires a conjugate factor. The factor r is 
interpreted as amplitude reduction due to tunnelling between 
edges. 
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Rigorously speaking, the second line is only heuristic, 
since the gauge vector potential A^ can not be smooth 
and single- valued everywhere on the torus (at r = 1) if 
JV^O; otherwise the 2-d integral in Eq. JBJ would van- 
ish due to Stokes' theorem. A rigorous treatment needs 
several patches to cover the torus, together with vector 
potentials pertaining to each patch. (For the details, see 
Re£M.) 



B. Proof of the Theorem 

To relate the bulk Chern number to edge states, we 
start with re-expressing the Chern number N. For any 
given held strength F xy on the torus at r = 1, we choose 
the gauge in which A x is single- valued and vanishes. 
Then we need two patches for A y to be single-valued 
in each. The simplest choice is 

AW(1,0 X ,9 V ) = / F xy (l,6 x ,9 v )d$' x ,0<e x <2ir, 

Jo 

A y 2 \l,9 x ,9 y ) = J F xy (l,9' x ,9 y )d9 x ,-n<9 x <it, 
A x (l,9 x ,9 y ) = 0. (7) 

Here A y ' and A y ' have discontinuity at 9 X = and 
9 X = it, respectively. In this way the torus is covered by 



two cylinders defined by 9 X ^ and 9 X ^ it. With this 
gauge fixing, the Chern number is expressed as 



N = 
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Defining 



2tt 



0(rA) = / d9 y A y 1 \r,9 x ,9 y ), 

Jo 
<£>(r,9 x ) = exp[i(f>(r,0 x )] , 

then Eq. JHJ can be written as 



2tt J 0+ 69 x 
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(9) 



2tt / r=1 
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Consequently, there is a 27riV jump in the 0-field, but the 
U(l) phase field $(r, 9 X ) is single- valued everywhere on 
the torus at r = 1. The Chern number is then equiva- 
lent to a winding number N of the 0-field on the loop 
r = 1,6*2, £ [0,27r). For more general gauge choices, the 
value of $ can change by a constant phase factor, but 
the relation between the Chern number and the winding 
number remains true as long as A x is kept to be single- 
valued. 

The physical meaning of the phase field <&(r, 9 X ) is sim- 
ply the Berry's phase that the system acquires due to an 
adiabatic evolution along the closed path with 9 y varying 
from to 27T and with r and 9 X fixed. Since the Hamil- 
tonian depends on r and 9 X through r e ±lTn:e:c only, all 
the parameters (0,0a:) corresponds to the same Hamilto- 
nian, which describes an open-boundary system. Conse- 
quently, the phase factor <&(r, 9 X ) is defined on a 2d plane, 
with (r, 9 X ) as the polar coordinates. 

Now the relation between Chern number and the gap- 
less edge states becomes clear. If N ^ 0, the non- 
vanishing winding number of (j>(r, 9 X ) on the unit cir- 
cle r = 1 requires vortex-like singularities in the disk 
bounded by the unit circle. Recall that the Berry phase 
$(r, 9 X ) is always well-defined by Eq. as long as the 
system 7J^ [c ^ y (' , ) 9 x ,9 y ) has a unique ground state and a 
gapped excitation spectrum for all 9 y E [0, 2tt). Thus we 
conclude that a singularity of the $-field can only occur 
at a point (r, 9 X ) where the system becomes either gapless 
or has ground state degeneracy for some 9 y . In summary, 
we have proved the following theorem: 

• Theorem 1 : Whenever the Chern number defined 
by Eq. © is non- vanishing, there must be a proper 



set of twisted boundary conditions {r,9 x ,9 y ) with 
r < 1, for which the Hamiltonian -ffr x r {T,9 x ,9y) 
has either ground state degeneracy or gapless exci- 
tations. 

As mentioned above, the boundary conditions for the 
system with r — 1 correspond to the usual definition of 
the twisted boundary conditions with twist phases 9 x ,9 yi 
while r — corresponds to defining a cylindrical system 
with open boundaries. The twisted boundary conditions 
with < r < 1 correspond to a (bent) cylindrical sys- 
tem with a certain amount of inter-edge tunnelling. The 
above theorem predicts only the existence of singularities 
of the Berry phase gauge field for systems with twisted 
boundary condition with r < 1, due to the appearance 
of gapless states, without providing detailed information 
about where the singularities are. It is possible that a 
system with a non- vanishing bulk Chern number has gap- 
less edge excitations in geometry with open-boundary (at 
r=0), like what happens in QH systems. But it is also 
possible that the gapless points shift away from r = 0, 
which means that the edge states are gapped in open- 
boundary systems but become gapless for some points 
(r 7 8 x ,9 y ) with r < 1. The latter case occurs for some 
QSH systems, as will be shown in next section. That both 
possibilities exist is the major new lesson we learn from 
Theorem 1, which improves our understanding of the re- 
lations between bulk topology and edge states. More 
precisely, we have learned that 

1. The bulk topological order, defined by first Chern 
number © and protected by the bulk gap, al- 
ways manifests itself by the existence of gapless 
edge states with a proper amount of inter-edge tun- 
nelling; 

2. The gaplessness of edge states in an open-boundary 
system is not protected solely by the bulk topology; 
its stability requires additional conditions. 



III. EXAMPLES: QUANTUM CHARGE AND 
SPIN HALL EFFECTS 



In this section, we will consider the usual integer quan- 
tum Hall effect and the recently proposed quantum spin 
Hall effect, as two simplest examples of topological insu- 
lators, and study the consequences of Theorem 1. For the 
quantum Hall effect, it is shown that the charge conser- 
vation plays an essential role in protecting gapless edge 
states in systems with open boundary. For the spin Hall 
effect the opposite occurs, i.e. the edge states in an 
open-boundary system may be gapped because spin is 
not conserved due to spin-orbit coupling. However, as a 
nontrivial consequence of the bulk topological order, the 
edge states can become gapless with increasing inter-edge 
tunnelling to a certain "critical" value. We will show this 
explicitly by numerical and analytic methods. 



A. Quantum Hall Effect: 
Role of Charge Conservation 

Compared with more general topological insulators de- 
fined by a nonzero bulk Chern number JSJl, the distinct 
feature of the quantum Hall effect is the charge conser- 
vation, i.e. the twisted boundary conditions are gauge 
equivalent to flux threading the torus and the Laughlin- 
Halperin gauge argument can apply. More generally, it 
is shown below that the edge states of an open-boundary 
system, -ffr x r {f = 0, x ,6 y ), are gapless in the thermo- 
dynamic limit if the following conservation conditions are 
satisfied: 



E 



c\v^a,H 



0. 



{n = x,y) 



[T x ,T y ] = 0, 



(11) 



where Hq is the untwisted Hamiltonian Q for a system 
with periodic boundary conditions. 

Under the conditions (flip, one can make a unitary 

transformation, U = exp <ij^ J^\ i x c\T x Ci >, to show 
that the twisted Hamiltonian J3J is gauge equivalent 
to the flux-threading Hamiltonian, -ffr x r (f,9 x ,9 y ) = 
UHt^t (r,9 x ,9 y )U^ . Namely H can be obtained by re- 
placing 



U 



%i j 6 



ir,. 



" \Jx tx ) 



V<y) 



(12) 



in the original Hamiltonian Hq. Here L is the size of 
the system along x direction. The second equation in 
Eq. Ulljl is necessary to maintain the conservation of 
T x ,T y in the twisted system. The boundary condition 
in y-direction is kept the same as in eq. |2")l. Conse- 
quently, the Hamiltonians H(r,9 x ,9 y ) and H{r,9 x ,9 y ) 
have identical energy spectrum. On the other hand, since 
A x — ► for L — ► oo, the energy spectrum of H should 
be insensitive to 9 X , except for the special case where the 
system has phase stiffness due to off-diagonal-long-range- 
order. In other words, if for a value ro ^ 0, the Hamilto- 
nian Hr x v (?"o, 9 x ,9 y ) were gapless in the thermodynamic 
limit, then Hp x r v (ro: 6' x i9y) should also be gapless for all 
9' x <G [0, 2tt). Thus if a singularity occurs for r ^ 0,9 X , 
then the points on the circle (with < 9 X < 2ir and the 
same ro) would all correspond to gapless systems, provid- 
ing a macroscopically large winding number. Therefore 
r must be zero. So we conclude that the gapless sin- 
gularity in parameter space for a system with a nonzero 
bulk Chern number N ^ in the thermodynamic limit 
can only occur at the point r = 0, corresponding to a 
system with open boundaries with gapless edge states. 

This argument, perhaps not very rigorous mathemati- 
cally, is enough to provide a better understanding of the 
role of charge conservation in the bulk-edge relation in 
the quantum Hall effect. In particular this implies that 
the stability of edge states in open-boundary systems is 
not protected by the bulk Chern number alone, if T^ or 



Ty are not generators of conserved quantities or if they 
do not commute with each other. This is exactly what 
happens in the case of the quantum spin Hall systems, 
which we will study in more details below. 



B. Quantum Spin Hall Effect: 
a Numerical Study 

The bulk topological number in the quantum spin Hall 
effect corresponds to choosing the internal generators in 
the twisted boundary condition (J2J to be 



(a) c=0, r=0 




s, 



1, 



(13) 



in the spin channel. Here S is a component of spin or 
pseudospin operators. There are two key differences be- 
tween the QSH and QH systems; namely (1) S may be 
non-conserved; (2) S is odd under time reversal (TR). In 
this section, we will focus on the consequences of point 
(1), and leave point (2) to next section. 

When S is not conserved, it is possible for the sin- 
gularity in the (r, 6 X ) plane to shift to some point with 
< r < 1. To see indeed this can happen in the twisted 
systems with edge tunnelling, below we will study the 
edge dynamics in a generalized version of the spin Hall 
system proposed in Ref. 1 -. The Hamiltonian for this sys- 
tem reads 



(b) o=0.05, r=0 
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in which T a ,a = 1,2,. .,5 are the generators of 5*0(5) 
Clifford algebra, and df 3 = -^ £ k d a (k)e ikr ' with d a (k) 
given by 

d\{k) = — vocsinfcy, 
(fo{k) = — v3csinfc x , 
dz{k) = — v3 sin k x sin k y , 

di{k) — vo (cos fez — cos k y ) , 

d$(k) = 2 — e s — cosk x — cos k y . (15) 

Since we are interested in the spin Hall insulator phase, 
direct hopping t is expected to be small. For simplic- 
ity, here and below we will take t — 0. In the case of 
inversion symmetry, d\ = di = and there is a con- 
served pseudo-spin operator: Ti2 = |Ti,r 2 ] /(2i). With 
the identification T^ = Ti2, T y = 1, the twisted boundary 
conditions J2J can be expressed as 



3(L,i/ 4 ),(l,l/0 



*(£,3/«),(l,j/«) re 



ie x r 12 



^d 



(xi,L), (xi,l) 



(16) 



Since I?i2 is conserved, the system is equivalent to a de- 
coupled bilayer quantum Hall system, and the bulk Chern 
number defined for the Hamiltonian Hr 12 ,i(r — 1, 6 X , 9 y ) 




(c) c=0.05, r=0. 18, 9 =n/2 




FIG. 3: The 1-d band structure for (a) c = and the open- 
boundary condition r — 0; (b) c = 0.05 and the open- 
boundary condition r — 0; (c) c = 0.05 and the twisted 
boundary conditions r = 0.18, 8 X — ±7r/2, 8 y = n. 



is A^ = 4 for < e s < 4, as shown in Ref^. Accord- 
ing to the arguments in the last subsection, there should 
be gapless edge states in the open-boundary system, as 
verified numerically in Fig[31(a). 

However, when di,d 2 7= 0, the edge states in the sys- 
tem with open boundary conditions will become gap- 
ful with a gap of the order \Jd\ + d|, since the r!,r 2 
terms will mix the states with opposite T^-eigenvalue. 
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FIG. 4: The distribution of (a) charge density; (b) pseudo- 
spin Fi2 density of the two edge states with wavevector k = 
0.16tt and with c = 0.05, r = 0.18, 9 X = tt/2. 



Also as long as |di|,|d2| <C \/d a d a ,a = 3,4,5, the 
bulk gap cannot be closed by turning on d\ , cfo ■ Conse- 
quently, the Berry phase gauge field remains well-defined 
on the torus with r = 1, and thus the Chern number 
for the Hamiltonian Hr 12 i(r — 1, x ,9y) remains un- 
changed. According to our Theorem 1, there must be 
some (r ^ 0,9 x ,8 y ) where the system become gapless. 
Since all the twist transformations only make changes in 
the boundary terms in the Hamiltonian, it is reasonable 
to expect the gapless excitations to live around the edges. 
The nonzero tunnelling r ^ between the edges actually 
mixes the original edge states to form new gapless ex- 
citations that are localized at the edge. To verify this, 
we have done a numerical calculation for the model with 
c = 0.05. We identified two gapless points in parameter 
space at r = 0.18, 6 X = ±ir/2,8 y = ir. The evolution of 
the band structure is shown in Fig. [3] At each gapless 
point there are two level crossings, which keeps the total 
winding number 4. The charge and pseudo-spin density 
distributions of the gapless states are shown in Fig. 0] 
confirming that it is localized symmetrically around the 
edge, and the pseudo-spin polarization is opposite on the 
two sides of the edge. 



C. Quantum Spin Hall Effect: 
a Perturbative Study 

It will be helpful to study how the gapless edge state 
arises from the original edge states through tunnelling 
in perturbation theory. Under the condition |di|, | c?2 1 ^ 
\/d a d a and r -C 1, all the effects induced by nonzero 
d\ , di and by twisted boundary conditions are only sig- 
nificant for the edge states, while leading to little change 
in the bulk states. Thus it is reasonable to describe the 
low energy dynamics in an effective 1-d edge Hamilto- 
nian. Let us first take the continuum limit of the gapless 
edge Hamiltonian for the system with d\ = d% = 0, r = 
as the starting point: 



m 



edge 



dkavk 



(4ls 



khsa ~ kRsa <j)kR Sa (17) 



in which a = ±1 corresponds to the eigenvalue of Ti2, 
and L, R the states at the left and right edges, s = 1, 2 
correspond to two crossings in Fig. [31 (a), and the wave- 
vector k is defined in reference to the corresponding cross- 
ing points. The time reversal transformation of the 8 



edge states is T^>, 



k,L{R),l,a 



T -l = 



a( P-k,L(R),2-a- 



When 



r / 0, we need to incorporate the mixing between dif- 
ferent edges, and when d±,d2 ^ the mixing between 
different a. The states with different s will never be 
mixed. So we take the perturbed effective Hamiltonian 
to be of the form 



H, 



edge 



+ r 



dkavk \4>t Lsa <t>kL B a ~ <t>{ Rsa 4>kRsc}j 
+ dk \ULs4>\ Ls] <i)kLsl + UR s <f>l m <f>kRsl + h.c.j 

E fdk{f kLsa e iea cl> kRsa +h.c.) 

s J 

£/ d **L M S*fc»i, (i,j = !,-••, 4) (18) 



in which $ fcs = (4>kLsi ,4>kLsi,4>kRs] , fikRsi)- The matrix 
elements ul(jq s are defined as 

u L{R)s = (k = 0, L(R)s T| cZiT 1 + d 2 r 2 \k = 0, L(R)s |) 
where the fe-dependence of the matrix elements is ig- 



nored. The matrix Mf, is defined as 



I vk 



'■Ls 



re 



M s 



° \ 



ULs 



vk re~ 1 ' 
re '" — vk u Rs 
\ re 10 u Rs vk J 



(19) 



Direct diagonalization of M gives the four eigenvalues: 



E sa p — a 



1 



v 2 k 2 + 



\ULs\ 2 + \uRs\ 2 



H 



\UR S \ 2 ~ \uLsf 



+ r 2 [\u Ls \ 2 + \u Rs \ 2 + 2\u Ls \\u Rs \cos29 s ) , (20) 



in which a, /3 = ±1 and 

9 S =9 + Arg(u Ls ) - Axg(u Rs ) . 

The equation E sa p = has two solutions: 

9* = or 7r, r = 1 



v\ 



ULsUR s \ 



(21) 



If r\ 7^ r<i, there are 4 gapless points (0, ri^), (7r, r\p) 
in the parameter space. In other words, the original sin- 
gularity at r = with winding number 4 is split to 4 sin- 
gularities each with unit winding number. However, the 
time reversal invariance will force r\ = ri-, since 4>kL(R)i 
are the T-partner of <p-kL(R)2- Therefore, the 4 singu- 
larities must merge to 2, each with winding number 2, 
which is consistent with the above numerical results. The 
schematic picture of <j) y (r, 9 X ) configurations for c = and 
c 7^ is shown in Fig. [SJ 

Concluding this section, we would like to emphasize 
that the existence of gapless edge states for a proper set 
of twisted boundary conditions with inter-edge tunnelling 
should not be viewed as a trivial result of fine tuning. For 
example, if one replaces the inter-edge tunnelling term in 
Eq. (QSJ by 



■E 



dk 



(d 
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kLsa e PkRsa 
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the Hamiltonian will be gapped for any (r, 9) as long as 
ul,Rs t^ and Im (ul s /ur s ) ^ 0, a direct consequence of 
the vanishing Chern number in the charge channel. 



IV. ROLE OF TIME REVERSAL INVARIANCE 

Time reversal invariance plays an important role in the 
topological description of the quantum spin Hall systems. 
For the non-interacting fermion case, a Z2 classification 
has been proposed in Ref^, namely the gaplessness of 
edge states in open-boundary systems is protected by 
Kramers degeneracy, if there are in total an odd num- 
ber pairs of edge states. However, the situation is very 
subtle when interactions are included in the many-body 
system. Below we will first show a consequences of time 
reversal invariance (Theorem 2) in a general many-body 
system in subsection A, and then discuss about the Z^- 
classification in subsection B, where new predictions will 
be made based on Theorem 1. 



FIG. 5: Schematic picture of the <j> y (r,8 x ) configuration for 
the spin Hall system when (a) c = 0; (b) c/0. The direction 
of the arrows stands for the phase cj> y . The crosses indicate 
the position of singularities and the circle is the unit circle 
r=l. 



A. Conditions for Winding Number Doubling 

For simplicity, we consider the case when T^ , T y each 
carries a simple representation of time reversal transfor- 
mation: 



Tr^T" 1 



V^„ 



Tfc = ±l. (22) 

If the original Hamiltonian Q is TR-invariant, then 
the TR property of the twisted Hamiltonian will be 

TH Tx v y {r,e x ,e v )T- x =H Ta: T v {r,-r) x 9 x ,-r} v ey) (23) 

Thus the gauge field has the symmetry property 

A M (r, 9 X , $ v ) = -1 (G(r, X ,0 V )\ 9 M \G(r, 9 x ,9 y )) 



= i (G(r, -rj x 9 x , -r) v d v )\dn \G(r, -r} x 6 x , -T] y 9 y )) 
= riptA^r, -rjxOx, -VyOy), 
=> F xy {r,6 x ,6y) = -rixVy^xyir, -r) x x ,-r) y 9 v ) . (24) 

Thus the Chern number N is non-vanishing only when 
VxVy = -l- I n the usual QHE case r] x — r) y = 1, 
which means that the Hall conductivity vanishes for TR- 
invariant systems. On the other hand, the spin Hall effect 
corresponds to rj x = — 1, rj v = 1, since T^ is a spin opera- 
tor and T y = 1 is taken to be the charge operator. Below 
we will focus on the case with r/ x = —l,rj y = 1. 

In general the first Chern number N can be either even 
or odd integer. However, we will show that N must be 
even when 



„ir x 7r 



-1 



(25) 



When T x is a component of the spin operator, this con- 
dition simply means that the spin is a half-odd integer. 

To prove this statement, it is convenient to introduce a 
deformed Hamiltonian by redefining the twisted bound- 
ary condition in the x-direction as 



U 



t i:j re ir * *e 10 



(ij) across L x 



(26) 



and leaving the twisted boundary condition in the y- 
direction unchanged as in Eq. JjJJ . In this way, we denote 
the auxiliary Hamiltonian as 



H Lr v (r,6 x ,a x ,6y). 



(27) 



The introduction of the charge twist phase e lQx changes 
the TR property of H' to 

TH^ ry (r,6 x ,a x ,9 y )T- 1 
= H^ Ty (r,9 x ,-a x ,-9 y ) . (28) 

On the other hand, noticing the condition 125|) . we get 

H v x r v (r,9 x ,a x ,9 y ) 
= H^ ry (r,9 x +TT,a x +TT,9 y ). (29) 

Combining Eqs. I|28|) and (|29|l and taking a x = n/2, we 
obtain 



TH^ Tv (r,9 x ,\,9 y )T-' 



(30) 



In the same way as for the Hamiltonian J2J, the 
Berry phase gauge field A' and the phase factor 

< t>'( r ^9 x ) — L A'(r,9 x ,9 y )d9 y can be defined for (the 
unique ground state of) H' . Then Eq. (|3l7|) leads to the 
symmetry property 



4>'(r,O x )=4'(r,0 x + n), 



(31) 



which means 4>' is symmetric under the central reflection 
(r,9) — > (r,9 + it). Consequently, the positions of the 



singularities are also symmetric under the same transfor- 
mation, thus the winding number must be even: 



N'(r) = 



d9 x 
* d<j>'(r,9 x ) 
d9 x 



d9 x 
d9 x 



even integer (32) 



if the Hamiltonian H^ r (r,9 x ,^,9 y ) is gapped for all 
9 x ,9 y € [0, 2n) at given r. 

As the last step, the no-phase-stiffness condition is in- 
troduced again, which means the energy spectrum of the 
system is insensitive to a x in the thermodynamic limit. 
Therefore, whenever H^ r (r,9 x ,^,8 y ) is gapped, the 
Hamiltonian J2J), 



Hr w r y (r,9 x ,9 y ) = H^ ry (r,9 x ,0,< 



y>- 



should also be gapped. This implies that the singularity 
distribution of cj)(r,9 x ) is the same as that of 4>'(r,9 x ). 
As a result, the Chern number of the spin Hall system 
Hr x r y (r, 9 X , 9 y ) must be an even integer, and we have 
finally proved the following theorem: 

• Theorem 2: If the conditions below are satisfied: 



iTTyT- 1 



r, 



1. TT X T- 1 = -IV 

2. exp(zr K 7r) = -1, 

3. the original Hamiltonian is TR-invariant, with 
charge conservation but without phase stiff- 
ness (i.e. without superconductivity), 

then the singularities in the phase field <fi(r, 9 X ) de- 
fined in Eq. @ distribute symmetrically with re- 
spect to the origin r — on the (r, 9 x )-pl&ne, and 
the Chern number for the unique ground state of 
ffr^j (r, 9 X , 9 y ) is an even integer. 

B. Comments on Zi Classification 

After studying the general consequence of time rever- 
sal invariance, now it is time to consider whether there is 
any additional topological classification besides the U(l) 
Chern number for TR-invariant systems. Although the 
singularities distribute symmetrically around the origin 
r = when the conditions in Theorem 2 are satisfied, it 
is still possible for all singularities to move away from the 
origin r = 0, which means the edge states in the open- 
boundary system is not protected to be gapless. However, 
it has been shown analytically 12,19 and numerically 15 
that gapless edge states along the open boundaries in a 
spin Hall insulator model proposed for graphene remains 
robust against disorder. In Ref>i£ it is further shown that 
the existence of gapless edge states and the non- vanishing 
bulk Chern number N = 2 in the spin channel are in 
one-to-one correspondence. The robustness of the gap- 
less edge states in an open-boundary system comes from 
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Kramers degeneracy, and applies only when there are an 
odd number of Kramers pairs on each edge, which corre- 
sponds to the Chern number N — An + 2, n e Z. In this 
way, a Z-i classification is well-defined for non-interacting 
fcrmion systems, which distinguishes the spin Hall insu- 
lators with Chern number N — An + 2 from those with 
N = An. 

To generalize this mechanism to an interacting system, 
it should be noticed that the condition for Kramers de- 
generacy, 



TTc ka T- l T- 1 



-Cko 



(33) 



is a single-particle property, not a condition on the many- 
body wave function. When interactions are taken into 
account, the gapless edge states emerging at some singu- 
lar point (r, X ) in parameter space may carry quantum 
numbers different from the constituent particle. Thus, 
the protection for Kramers degeneracy can survive only 
if the condition \\.\'M is still true for the edge modes. Such 
a protection can be understood both by edge dynamics 
and by bulk topology. Below we will discuss about Zi 
classification from the bulk topology as a consequence of 
both Kramers degeneracy and Theorem 2. 

To proceed, we introduce a conjecture that is known 
to be true for non-interacting systems^, but has not been 
proven for a general many-body system: 

• Conjecture: If there is a singularity of the <fi(r, X ) 
field with winding number N at some point (r, X ), 
the corresponding Hamiltonian H-p^r (r, X , y ) has 
2JV branches of gapless edge states. 

Below we will take this conjecture as the starting point 
and leave its proof to future work. In accordance with 
this conjecture, a singularity with winding number N = 1 
corresponds to 2 branches of edge modes. If the edge 
modes carry T 2 = — 1, the two branches must be T- 
conjugate to each other. Therefore there must be a 
Kramers degeneracy between them that cannot be bro- 
ken by any perturbation without breaking the time rever- 
sal symmetry. Recall that the twist phase factor e lTx " is 
T-even when TT x T~ l = —T x , the local variation of pa- 
rameter (r, X ) around the singular point can be consid- 
ered as local perturbation respecting T-symmetry, which 
thus must maintain the system gapless. In other words, 
the system cannot have point-like singularity with wind- 
ing number N = 1 anywhere in the parameter space. A 
similar argument shows that all the point-like singular- 
ity in the (r,0 x )-pl&ne with N odd is forbidden. Con- 
sequently, all the vortex- like singularities of the 4>( r ,0 x ) 
field has an even winding number. Taking into account 
the central reflection symmetry of the singular points as 
a corollary of Theorem 2, we know that all the singulari- 
ties away from r — provide a winding number An. Thus 
at least 2 pairs of gapless edge states must exist for an 
open-boundary system when the bulk Chern number is 
N = 4n + 2, (ne Z). 

In summary, the Zi classification is applicable when 
the system satisfies all the conditions of Theorem 2 and 



also has all its low- energy quasiparticles odd under T 2 - 
operation. It should be clarified that for an interact- 
ing system, the condition (|25|l does not imply T 2 = — 1 
for low- lying modes. Since the properties of low-energy 
quasiparticles may change without closing the bulk gap, 
there seems no reason to believe the Zi "order" is pro- 
tected by the bulk topology alone; it rather depends on 
the detail of interactions. 

From the discussion above, new predictions can be ob- 
tained for the N = 1 system. In Refi22*2I, the edge dy- 
namics for some N = 1 systems has been studied, which 
shows that the edge states in the open-boundary system 
can be gapped with proper interactions. After such a gap 
is opened, there are two possibilities for this system: 

1. As shown in Ref^S, the ground state of an 
open-boundary system may break time reversal 
symmetry spontaneously, which leads to ground 
state degeneracy. In this case, the point r = 
is still singular. If the Chern number N = 2 
is still carried by this singularity, then no other 
singularity necessarily exists. 



2. If the singularity at r = does not carry any Chern 
number, or if the TR-symmetry is recovered, then 
the bulk topology requires the gapless edge states 
to exist with proper inter-edge tunnelling strength. 
What's more, from the discussion in this subsection 
we know that as long as the gapless edge states shift 
to r 7^ 0, they must be even under T 2 operation. 
Physically, this prediction means that before the 
gap-opening transition, the edge modes are T 2 -odd 
quasi-particles, which carry similar quantum num- 
bers as free fermions; immediately after this tran- 
sition, the edge modes will become T 2 -even, which 
may correspond to some kind of spin and/or charge 
density waves. 



V. SUMMARY AND DISCUSSIONS 

In conclusion, in this paper we have established, for 2d 
insulators, a general connection between bulk topological 
order and edge dynamics illustrated by Theorem 1. The 
key idea behind demonstrating this bulk-edge relation is 
to bring in an additional parameter that describes the 
strength of edge tunnelling besides the twisted bound- 
ary phases introduced in Refii This results in a three 
dimensional parameter space and implements a continu- 
ous interpolation between a cylindrical system with open 
boundaries and a toroidal system. In this framework, it is 
clarified that generically the bulk topological order for an 
insulator in torus geometry is not necessarily associated 
with gapless edge states in cylinder geometry. On the 
one hand, the the emergence of gapless edge states known 
for QH systems is shown to be closely related to charge 
conservation in addition to topological reasons. On the 
other hand, a QSH system with non-trivial bulk Chern 
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number and gapped edge states is predicted to have gap- 
less edge states when inter-edge tunnelling with proper 
strength is turned on. Finally, a general consequence of 
time reversal invariance is deduced in Theorem 2. The 
Z^ classification is shown to be dependent, besides bulk 
topology, also on the T 2 (T being time reversal) quan- 
tum number of relevant low-lying modes, which is not 
topological in nature and can be changed by interactions 
without changing the bulk topological order. 

The proof of Theorem 1 can in principle be general- 
ized to topological insulators in higher dimensions, for 
which the topological order is related to non-Abelian 
Berry phase gauge field and higher Chern numbers. For 
example, in the four dimensional quantum Hall effect 24 , 
the second Chern number can be defined when the Berry 
phase gauge field is SU(N) and the parameter space is 
4-dimensional as labeled by {0 x ,8 y ,9 z ,6 w ). In the same 
way as in Sec. II, an amplitude r can be introduced 
for any boundary and a non-Abelian Berry phase field 
U(r, 0,j,) £ SU(N) can be defined by a straightforward 
generalization of Eq. JJjJ to a Wilson loop operator: 



$(r,6^) =Pexp 



i / A w (r,9 x ,8y,8 z ,9 w )d8 u 



in which \x = x,y, z and P stands for path ordering. The 
second Chern number is then identified as the winding 
number of the U field on the 3-torus T 3 at r = 1. By the 
same topological reasoning, there must be singularities 
of the U field for some (r, 6^) with r < 1 when the bulk 
Chern number for the system with r = 1 is non-trivial, 
which corresponds to gapless edge states in higher dimen- 
sional sense. In this way, Theorem 1 provides a universal 
modus operandi to understand the bulk-edge relationship 
in topological insulators. 



(34) 
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